Organic Molecules on Wide-Gap Insulators: Electronic Excitations 
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The electronic excitation of a conjugated molecule- insulator interface, as exemplified by the ad- 
sorption of benzoic acid and its phenolic derivative on NaCl(OOl) surface, is addressed by many-body 
Green's function methods. By solving the two-particle Bethe-Salpeter equation on top of the GW 
quasiparticle energies, it turns out that instead of the intramolecular ty — tv* transition of the adsor- 
bate, the lowest singlet excited state of the adsorbate system is essentially assigned to the transition 
from the surface valence band maximum to the 7r* state of the molecule, as a charge transfer exci- 
tonic effect. An accurate description of this lowest electronic excitations confined at the interface 
requires the knowledge of a full excitonic Hamiltonian due to the sizable electron-hole exchange 
interaction. 

PACS numbers: 71.35.Cc, 73.20.-r, 78.20.-e 



A revived interest in the organic molecule and wide- 
gap insulator interfaces has been witnessed in the last 
decade. Their promising applications in molecular elec- 
tronics [1] are benefited from the merit that the electronic 
characteristics of the adsorbed molecule are barely per- 
turbed by the underlying insulating surfaces [2[, which 
enables the use of wide-gap insulator as a supporting sub- 
strate for a wide variety of technical applications. The 
organic molecule-insulator interface is also relevant to the 
separation process seen in mining industry, namely vari- 
ous mineral species, usually rocksalt (NaCl) and sylvite 
(KC1) can be distinguished electrostatically when mixing 
with certain acids, e.g. salicylic acids (SA). It is specu- 
lated that in the combined molecule-insulator system, the 
presence of the highest occupied (HOMO) and the lowest 
unoccupied molecular orbitals (LUMO) can enhance the 
charge transfer by electronic excitations by several orders 
of magnitude with respect to the bare insulators, because 
the molecular energy gap is usually much smaller than 
the band gap of the wide-gap insulators [3|, |4| . When the 
adsorbate covered insulating surfaces are brought into 
contact, those excited electrons, as mobile electrons, can 
hop in-between the surfaces until the thermal equilibrium 
is reached and the Fermi levels are aligned across the in- 
terface. This gives rise to a contact voltage due to the 
charge transfer, which tentatively explains the origin of 
the electrostatic separation. 

Previous Kohn-Sham density functional theory (KS- 
DFT) calculations 0, ^ have identified that the ad- 
sorption of benzoic acid (BA) and its phenolic deriva- 
tives (e.g. SA) on alkali chloride surfaces indeed yields 
a much smaller effective KS gap than that of the bare 
surface. However, the KS gap, defined as the difference 
between the eigenvalues of the lowest unoccupied and 
highest occupied KS states, is not a realistic measure 
of excitation properties. One of the basic issues with 
KS gap is that it is always smaller than the fundamen- 
tal gap by a derivative discontinuity of the exchange- 
correlation energy functional [7] even if the discontinu- 



ity is present in the functional [8|, and this error could 
be sizable for wide-gap insulators [9]. This gap problem 
can be alleviated either by hybrid functionals incorpo- 
rating nonlocal Fock exchange within the generalized KS 
scheme [10| , or ultimately from the quasiparticle energies 
in the GW approximation of the self-energy [ll| . On the 
other hand, in context of the electronic excitations at 
the molecule-insulator interface, the electron- hole (e-h) 
interaction is not directly accessible from the fundamen- 
tal gap within a single-particle picture, but only available 
in time-dependent DFT [l2| and the two-particle many- 
body perturbation theory (MBPT), namely the Bethe- 
Salpeter equation (BSE) approach [13j. The excitonic ef- 
fect has been found prominent for some organic molecules 
and wide-gap insulator surfaces [17j in terms of 
MBPT. Nevertheless, little is known about the electronic 
excitations when organic molecules are adsorbed on in- 
sulating surfaces. In this Letter, we demonstrate that, 
based on Go Wo and BSE calculations of BA and SA ad- 
sorbed on NaCl(OOl) surface, the lowest singlet excited 
state for a 7r-conjugated molecule-insulator interface in- 
volves the contributions from both the surface and the 
adsorbate molecule, instead of the intramolecular tt — 7r* 
transition. This is not readily expected provided that the 
coupling of the molecule and the surface is weak. 

The GoWo and BSE calculations were performed on 
top of the KS-DFT eigenvalues and eigenstates within 
local density approximation (LDA). The ground-state ge- 
ometries were also optimized using DFT-LDA. We em- 
ployed plane-wave basis sets with a cutoff energy of 70 
Ry, with Troullier-Martins norm-conserving pseudopo- 
tentials as implemented in abinit [18]. The LDA energy 
levels were corrected to the first order in the subsequent 
GoWo calculations, where the screened Coulomb inter- 
action W is expressed in random phase approximation 
(RPA) (W — s~ x v\ with e and v being the dielectric 
function and the bare Coulomb interaction. The inverse 
dielectric function e~ x was evaluated in the plasmon-pole 
approximation [19j. Spurious Coulomb interactions aris- 
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ing from the periodic images were truncated by a box-like 
cutoff [20]. The excitonic effect is finally taken into ac- 
count by solving the BSE using the quasiparticle energies 
from the Go Wo calculations, which is equivalent to find- 
ing the eigenvalues of an effective two-particle excitonic 
Hamiltonian H exc as 



-[H c » 1 ]* 



-[H res ] 



(1) 



In terms of the one-particle occupied and empty quasi- 
particle states (v,c) and quasiparticle energies £^ P , the 
Hermitian resonant term i^ res , containing the transitions 
from the occupied to the empty states (v — > c) of posi- 
tive energies, can be expressed as (s^ p — £® P )S VV '5 CC ' + 
(K d + 2if x )( vc )(y c /), where K d and K x refer to the 
(screened) direct e-h Coulomb interaction and the bare 
e-h exchange interaction, respectively. The direct in- 
teraction is calculated with the static screening in this 
work. Analogously, the anti-resonant part — [i^ res ]* con- 
tains the c — >• v transitions of negative energies. The off- 
diagonal term H cpl = (K d + 2K X )( VC ^ C , V ^ describes the 
coupling between the resonant and anti-resonant transi- 
tions. Within the Tamm-Dancoff approximation (TDA), 
this coupling part is often neglected for optical absorp- 
tion calculations since it is usually much smaller than 
the resonant part [2l[. The TDA significant reduces the 
computational cost and becomes the routine for the BSE 
calculations for the excitations in various systems. How- 
ever, it is recently found that the overestimation of the 
excitation energy by TDA can be significant for organic 
molecules and nanostructures with confined excitations 
15l . fl6l . [22 1 . With an efficient iterative approach to the 
non- Hermitian H exc [16], it is thus timely and necessary 
to go beyond the TDA in the present work. The Go Wo 
and BSE calculations are carried out using the yambo 
code 23(. 

We first consider the gas-phase molecules. BA is the 
simplest aromatic carboxylic acid, with one carboxylic 
group attached to an aromatic ring. BA can be func- 
tionalized to SA when one phenolic group is attached 
to the ortho position. We show in Fig. [T] the molecu- 
lar orbital energies calculated by LDA and G0W0-LDA. 
The molecules are confined in a box with the dimen- 
sion of 30.0 x 30.0 x 24.0 Bohr 3 . The k-point is sam- 
pled at the T point. Over 450 empty states are included 
for the response function (with a 4 Ry cutoff) and the 
Green's function for the self-energy (with a 68 Ry cut- 
off for the exchange self-energy). It is not surprising 
that LDA distorts the occupied molecular orbitals when 
compared to the experimental ultraviolet photoemission 
spectra (UPS) [4]. Notably, LDA gives a different or- 
der of the three highest occupied orbitals of BA with re- 
spect to the Go Wo result. The Go Wo predicts that both 
the HOMO and HOMO-1 are delocalized ir states of the 
cyclic benzene ring, whereas it turns out that the LDA 
HOMO is the non-bonding state localized around the car- 
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FIG. 1. Molecular orbital energies of BA and SA calculated 
with LDA and G0W0-LDA, compared with experimental UPS 
data taken from Ref. Q. The highest occupied orbitals are 
aligned to zero eV. The wavefunctions of the frontier orbitals 
by the PW1PW hybrid functional are presented. 



boxylic group, which corresponds to the HOMO-2 in the 
G0W0 calculation. This is interpreted as a consequence 
of the spurious self-interaction in the local and semilocal 
approximations of KS-DFT, which tends to shift the lo- 
calized electronic states to higher energies |8[ . The effect 
of self-interaction is also observed for SA, where the LDA 
HOMO-2 is blue-shifted towards the HOMO-1 by about 
0.6 eV compared to the quasiparticle energies. We note 
that the self-interaction error can be effectively lifted by 
the nonlocal exact-exchange in hybrid functionals, and 
our PW1PW hybrid functional calculations [5] are in line 
with the Go Wo quasiparticle energies for the occupied or- 
bitals. Being a dynamic theory, the G0W0 results are in 
quantitative agreement with the UPS experiment [4[ . In 
particular, the G0W0 HOMO-3 levels for both molecules 
align well with the UPS peaks. 

Moving to the unoccupied orbitals, we find that the 
quasiparticle energy corrections are rather substantial. 
An interesting observation is that the quasiparticle cor- 
rection to the LUMO+2 is about 2 eV smaller than 
that to the LUMO and LUMO+1 in the present calcu- 
lations. Instead of the 7r* anti-bonding character of the 
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TABLE I. The lowest singlet excitation energies (in eV) of 
BA and SA calculated within BSE, in comparison with exper- 
imental values. The exciton binding energy Eb is calculated 
as the difference between the G W HOMO-LUMO gap (E g ) 
and BSE excitation energy. 
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TABLE II. Calculated tt - tt* gaps (in eV) of BA and SA on 
NaCl(OOl) by LDA and Go Wo -LDA. The isolated molecule 
refers to the molecule detached from the surface while its ge- 
ometry is kept fixed as that of the adsorbate. The molecular 
energy gap change upon adsorption is denoted by AE S . 
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FIG. 2. Optical absorption spectra of BA and SA obtained 
by full BSE calculations. The spectra are broadened by an 
artificial Lorentzian of 0.1 eV. 



LUMO and LUMO+1, the wavefunction of the LUMO+2 
is largely delocalized away from the molecule, a sign of a 
Rydberg orbital |24[. The Rydberg orbital is in general 
not correctly described by LDA because of the unrealis- 
tic exponential decay of the potential in the long-range. 
In terms of the HOMO-LUMO gap, we see in Fig. CD that 
the GoWo quasiparticle energy gap is more than twice as 
large as the LDA gap. 

Referring to the lowest absorption peaks from available 
26| . one notices that the GoWo 



UV/ visible spectra [25 



gaps, as the onsets in the GW-RPA independent-particle 
spectra, are also much larger than the experimental ex- 
citation energies (cf. Table H]), implying that the bound 
e-h interaction must play an important role in the optical 
absorption if the GoWo gaps are reasonable. When the 
e-h interaction is taken into account in the BSE, excel- 
lent agreement with experiment can be achieved in terms 
of the lowest excitation energy. In the BSE calculations 
for the molecules, the cutoff energies for the direct (K d ) 
and exchange (K x ) e-h interactions are 2 Ry and 68 Ry, 
respectively. We include e-h pairs involving the transi- 
tions in an energy window of 33.2 eV in the BSE kernel, 
as shown in Table HI It is seen in Table U that the devia- 
tion between the BSE excitation energy and experimental 
optical peak is within 0.12 eV. The coupling part in the 
excitonic Hamiltonian has a small effect on the lowest 
singlet excited state for the individual molecules, as the 
peak shifts by merely 0.1 - 0.2 eV using the TDA. The 
optical absorption in Fig. [2] is given by the imaginary 
part of the macroscopic dielectric function £m(q = 0,cj), 



where q is the momentum transfer. The lowest singlet 
excited state for both BA and SA largely correspond to 
the G Wo HOMO LUMO transition, i.e. the tt tt* 
transition. The polarization direction dependence of the 
lowest excitation peak can be understood by the symme- 
tries of the HOMO and LUMO wavefunctions illustrated 
in Fig. [TJ The e-h pair is very strongly bound, as the 
molecular exciton binding energies shown in Table [I] are 
about half of the fundamental gaps. Due to the nature of 
the 7r — » 7r* transition, the exciton is spatially localized 
at the aromatic benzene ring. The strong excitonic effect 
is an important feature of the benzoic and its phenolic 
derivatives. 

After examining the individual molecules, we now pro- 
ceed with the adsorption on the NaCl(OOl) surface. The 
molecules are placed on a four-layer (2x2) NaCl(OOl) su- 
percell, with a vacuum thickness of at least 15 A. The 
LDA bulk lattice parameter (a = 5.64 A) is used. We 
use a T centered 3x3 Monkhorst-Pack k-point mesh (5 
k-point in the irreducible Brillouin zone). The adsorp- 
tion geometry optimized by the LDA functional is justi- 
fied since the long-range van der Waals (vdW) interaction 
has a negligible influence on the adsorption configuration 
[6|. The molecules are weakly bound to the surface, with 
a binding energy of about 0.4 eV from the ionic contribu- 
tion and 0.5 eV from the vdW interaction [6]. While the 
unit cell is small, the lateral intermolecular interaction is 
still limited, as is manifested by the small energy level dis- 
persion of the HOMO and LUMO with respect to the k 
points. For GoWo calculations, we applied a total number 
of 960 bands in evaluation of the self-energy, and a cutoff 
of 60 Ry and 4 Ry for the exchange self-energy and the 
response function, respectively. Since the lowest excited 
state of the molecules involves the tt — » 7r* transition, the 
evolution of the tt — 7r* energy gap upon adsorption on 
NaCl(OOl) was investigated and is presented in Table HI1 
for BA and SA. Note that in LDA the tt — 7r* gap refers to 
the (HOMO-l)-LUMO instead of the HOMO-LUMO gap 
for BA. Analogous to the gas-phase molecules, the GoWo 
energy gaps of the adsorbates are significantly larger than 
the LDA ones. In particular, we find that LDA predicts 
that the HOMO of BA is below the surface valence band 
maximum (VBM), which is clearly at variance with our 
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FIG. 3. Absorption spectra of the BA and SA molecules ad- 
sorbed on NaCl(OOl) calculated within the Tamm-DancofT 
approximation (TDA) and full excitonic Hamiltonian. The 
perturbative polarization direction with respect to the sur- 
face is denoted by 0. The spectra are broadened by a 0.1 eV 
Lorentzian. 



earlier hybrid functional results |5[ as well as with the 
present Go Wo calculation. As the non-bonding state lo- 
calized at the carboxylic group is in resonance with the 
surface VBM [5|, the misalignment of the energy lev- 
els in the LDA calculation is clearly associated with the 
self-interaction error, which places the non-bonding state 
too high in energy. It is further revealed in Table [XT] 
that, while LDA yields almost identical energy gaps for 
the isolated and the adsorbed molecules, the quasipar- 
ticle energy gaps are reduced by about 1 eV when the 
molecules are in contact with the surface according to the 
Go Wo calculations. The renormalization of the molecu- 
lar energy gap in Go Wo is essentially an effect of the 
polarization in the surface due to the added or removed 
electron in the adsorbate, which in turn acts back on 
the adsorbed molecule and shifts its electron affinity and 
ionization energy level positions [27j. For example, the 
ionization energy level (HOMO) moves up by 0.36 eV 
whereas the electron affinity level (LUMO) moves down 
by 0.71 eV when a SA molecule is adsorbed on NaCl(OOl). 
The asymmetric shifts suggests that the renormalization 
in the present case does not strictly follow the classical 
image potential theory, and polarizations between the 
molecules are likely to take place with a full monolayer 
coverage. Decomposing the self-energy into the bare ex- 
change (E x ) and the correlation (S c ) components, we 
find that the shifts of the energy levels arise solely from 
the change in the E c . The polarization effect, thus being 
a dynamic correlation effect, cannot be captured by the 
KS eigenvalues in standard DFT approximations. Our 
results show that the effect can be pronounced even on 
wide-gap insulators with small dielectric constant, such 
as NaCl. 



As the final step, the excitonic effect at the molecule- 
surface interface is addressed by BSE. We include e-h 
pairs with energies up to 33 eV in the two-particle re- 
sponse function. The cutoff energies for the K d and K x 
were chosen to be 3 Ry and 20 Ry, respectively. The 
lowest excitation energies are converged within 0.2 eV. 
In Fig. [3] the optical spectra calculated within the TDA 
and the full BSE calculation are presented. We first dis- 
cuss the TDA results. The first excitation peak appears 
at around 4.1 eV for both the BA and SA adsorbate sys- 
tems. This peak is assigned to the tt — >• 7r* transition 
of the molecule according to a full diagonalization of the 
TDA Hamiltonian. Going from the parallel perturbing 
field (6 = 0°) to the perpendicular polarization (0 = 90°), 
the intensity of the 7r — > 7r* peak does not change much 
because the molecules are tilted on the surface. Moving 
to higher photon energies, we see that the second low- 
est peak appears at 4.6-4.7 eV for both molecules, which 
primarily comprises the excitation from the surface VBM 
to the molecular 7r* orbital, i.e. the CI 3p — > 7r* tran- 
sition. This peak is thus clearly attributed to a charge 
transfer exciton, which in principle arises from the direct 
overlap between the surface CI states and the delocalized 
7r* orbital of the adsorbate. The observed charge trans- 
fer exciton is rather surprising, because in the ground 
state the charge transfer is very small from the surface 
to the molecule {^j. At 6.0-6.1 eV, the strongest absorp- 
tion peak appearing in the parallel polarization stems 
from the NaCl(OOl) surface exciton, which lies about 1 
eV below the bulk exciton due to the image states outside 
the surface [17] . The peak related to the surface exciton 
is suppressed for the transverse field. 

When the coupling part is switched on in the full BSE 
calculation, the valence band 3p — » 7r* peak experiences a 
sizable red-shift by 1.3-2.0 eV. Remarkably, not only the 
excitation energy shifts, but the intensity of the 3p — » 7r* 
transition also grows much stronger when = 90°. Such 
a large perpendicular transition dipole moment further 
evidences that the exciton is spatially distributed along 
surface normal. The much stronger oscillator strength of 
this peak for the adsorbed SA molecule is benefited from 
its adsorption configuration, as the SA molecule is more 
parallel and closer to the surface than the B A [6| , giving 
rise to a stronger overlap between the VBM and the 7r* 
orbital. On the other hand, the lowest singlet molecular 
excited state and the surface excited state are not much 
influenced by the coupling part. As a result, the lowest 
excited state for the adsorbate system is now replaced by 
the 3p —> 7r* transition, with excitation energies of 3.3 
eV and 2.7 eV for the BA and SA adsorbate systems, 
respectively. Our results shows that the presence of the 
charge transfer exciton considerably lowers the effective 
excited energy of the molecule-insulator system, which 
may serve as the key ingredient in the contact charging 
model. 

The striking failure of TDA reveals that the coupling 
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part becomes comparable to the resonant part in the ex- 
citonic Hamiltonian for the 3p —> 7r* transition. In fact, 
we find that even without the direct interaction K d , the 
coupling term with the exchange interaction K x alone 
already shifts the lowest excitation energy from 4.6 eV in 
TDA to 3.0 eV in the case of SA. Therefore, the 3p tt* 
exciton has a prominent exchange interaction between 
the electrons and holes, making TDA insufficient in de- 
scribing this excited state. 

To conclude, we find that even when a conjugated 
organic molecule is weakly adsorbed on an wide-gap 
insulating surface, the presence of the exciton at the 
molecule-insulator interface can provide a much more 
effective channel for electronic excitations than the in- 
tramolecular transitions. The formation of the lowest 
excited state involving transitions from the surface to 
the molecule is driven by the sizable e-h exchange inter- 
actions. 

This work was supported by K+S AG. We acknowledge 
generous computation time allocated at HLRN. 
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